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ABSTRACT 

Theoretical studies and current observations of the high-redshift intergalactic medium (IGM) indicate 
that at least two cosmic transitions occur by the time the universe reaches gas metallicities of about 
10 _3 Z Q . These are the cosmological reionization of the IGM, and the transition from a primordial to 
present-day mode of star formation. We quantify this relation through new calculations of the ionizing 
radiation produced in association with the elements carbon, oxygen and silicon observed in Galactic 
metal-poor halo stars, which are likely second-generation objects formed in the wake of primordial 
supernovae. We demonstrate that sufficient ionizing photons per baryon are created by enrichment 
levels of [Fe/H] ~ -3 in the environment of metal-poor halo stars to provide the optical depth in the 
cosmic microwave background of ~ 0.1 detected by WMAP. We show, on a star by star basis, that a 
genuine cosmic milestone in IGM ionization and star formation mode occurred at metallicities of 10 -4 - 
10~ 3 Zq in these halo stars. This provides an important link in the chain of evidence for metal-free first 
stars having dominated the process of reionization by z ~ 6. We conclude that many of the Fe-poor halo 
stars formed close to the end of or soon after cosmological reionization, making them the ideal probe 
of the physical conditions under which the transition from first- to second-generation star formation 
happened in primordial galaxies. 

Subject headings: cosmology: theory — galaxies: high-redshift — nucleosynthesis - stars:abundances 
— stars: Population II 



1. INTRODUCTION 

Theoretical studies of the early universe predict the 
existence of first-generation stars which lack metals and 
consequently have unique physical structures and proper- 
ties (Bromm et al. 2001b; Tumlinson et al. 2003; Schaerer 
2002). Such objects have not, however, been detected to 
date in either the high-redshift or local universe. There- 
fore, indirect inferences from modelling the impact of pri- 
mordial stars and supernovae (SNe) on their environment 
through radiative and chemical feedback remain impor- 
tant. The reionization of the intergalactic medium (IGM) 
and the widespread metal enrichment of high-z galax- 
ies, QSOs and the IGM provide strong observational con- 
straints for such theoretical studies. A particularly effec- 
tive tool in constraining the masses and properties of the 
first stars has been to combine their ionizing photon and 
nuclcosynthctic signatures (Venkatesan & Truran 2003; 
Tumlinson et al. 2004). 

When we examine the current data and the results of 
theoretical work, an interesting coincidence becomes ap- 
parent: that three important cosmic transitions occur at 
gas metallicities of about 1O _5 -1O _3 Z , 0. To begin with, 
the universe first generates sufficient ionizing photons for 
hydrogen reionization when the IGM reaches these metal- 
licities (Miralda-Escude & Rees 1997; Venkatesan & Tru- 
ran 2003). Second, theoretical studies of the cooling and 
fragmentation of gas in early galaxies reveal that the pri- 
mordial stellar initial mass function (IMF) evolves from its 
metal-free to present-day form at gas transition metallici- 
ties of Z tt ~ 10- 5 -10~ 3 Z Q (Bromm et al. 2001a; Schnei- 



der et al. 2002; Omukai et al. 2005; Santoro & Shull 
2006). This transition results from the physics of gas cool- 
ing in early galaxies, when the dominant cooling chan- 
nel switches from molecular hydrogen to metal lines or 
dust grains, with the latter pathway dominating in sys- 
tems with Z < 1Q- 4 Zq (Omukai et al. 2005). Third, there 
are indications that the relative abundances of CNO and 
alpha-elements in metal-poor stars in the Galactic halo 
have increased scatter and altered ratios when the iron 
content of these stars falls below ~ 10~ 3 of solar values 
(Barklem et al. 2005; Cayrel et al. 2004; Beers & Christlieb 
2005, and references therein). These subsolar mass ob- 
jects are thought to be second-generation stars that formed 
from metal-poor fragmenting gas in the shells of the very 
first SNe in early galaxies (Salvaterra et al. 2004). Such 
stars contain the fossil imprint of nucleosynthesis from the 
preceding generations of metal- free stars, reflecting the 
"prompt inventory" of metals (Qian & Wasserburg 2002; 
Wada & Venkatesan 2003). 

When these three points are taken together, we see that 
several cosmological milestones occured by metallicities of 
~ 10~ 3 Zq, a relation that must result from the strong 
feedback between the buildup of ionizing radiation and 
the formation of second-generation stars. In this Letter, 
we quantitatively demonstrate a novel connection in the 
environment of extremely metal-poor (EMP) halo stars 
between two of these relationships: that the stellar popu- 
lation that was important for early nucleosynthesis is also 
responsible for cosmological reionization. We define EMP 
stars as those with [Fe/H] < -3 (Beers & Christlieb 2005). 
We examine the significance of the second transition in 
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relation to the other two in a subsequent work, where we 
plan to use numerical simulations to investigate the spa- 
tially varying feedback between gas cooling and the prop- 
agation of ionizing radiation in the first galaxies. 

2. STELLAR ABUNDANCES AND IONIZING EFFICIENCY 

We begin by showing the trends of element abundances 
in EMP stars from current observations. We follow our ap- 
proach in Venkatesan et al. (2006) and take the data from 
Cayrel et al. (2004), and Tumlinson et al. (2004) and refer- 
ences therein. For the two most iron-poor stars, HE-0107- 
5240 and HE-1327-2326, we use the published abundances 
from Christlieb et al. (2004), Aoki et al. (2004), Frebel 
et al. (2005) and Frebel et al. (2006). The last of these 
papers has recently revised the CNO abundances of HE- 
1327-2326 when 3D NLTE corrections are taken into con- 
sideration. We focus on the elements C, O, Si and Fe, as 
these are amongst the dominant metals created in primor- 
dial SNe and are the most important coolants in metal-free 
gas (Bromm & Loeb 2003; Santoro & Shull 2006). 

The trends of [C/H], [O/H], [Si/H] and [C/O] as a 
function of iron abundance [Fe/H] are shown in Figure 
1, where for an element X, we define [X/H] as the ratio 
of the measured column densities of X/H to its solar ra- 
tio. We display those points with [Fe/H] < -2 in order to 
span the predicted Z tr in individual elements (Bromm & 
Loeb 2003; Santoro & Shull 2006), as we will see below. 
We assume that EMP stars with [Fe/H] < -3 have been 
enriched only by metal-free stars, consistent with the esti- 
mate of Z tr from Santoro & Shull (2006) and with recent 
calculations using merger tree models in which stars of 
metallicities < 1Q~^Zq reflect enrichment from only one 
to a few primordial SNe (Tumlinson 2006). 
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Fig. 1. — The current data from Galactic metal-poor halo stars 
on the abundances of C, O and Si relative to H, and of [C/O], as a 
function of [Fe/H]. 

Figure 1 reveals some interesting trends. First, for those 
stars that have both abundances measured, [C/O] shows 
no clear trend with [Fe/H]; in fact, it appears to have its 
maximum value at the lowest [Fe/H]. Such element abun- 
dance ratios are often used as a cosmic clock to track the 
creation of elements from stars of differing lifetimes, but 
such a relation may not be easily derived in EMP stars, a 
point which we return to in the discussion. Second, car- 



bon and silicon roughly track Fe in enrichment relative 
to H, whereas oxygen does not. Again, the exceptions to 
this trend are the most Fe-poor stars. Note that some of 
the stellar data shown here come from surveys that se- 
lect against C-rich EMP stars (Cayrel et al. 2004), owing 
to the complicated analyses involved in such cases. The 
more recent survey by Barklcm et al. (2005) does not have 
this bias, although these authors focus mostly on r-process 
element abundances in EMP stars. 

A direct translation of the observed abundances in EMP 
halo stars to yields from primordial SNe and constraints 
on the first-stars IMF would be ideal, but is complicated 
by several factors. Mass transfer from a companion may 
alter original abundance ratios, but is difficult to model 
owing to the unknown binary fraction of the first stellar 
generations, and the role of AGB phase dredge- ups and 
mass loss for metal-free stars. We have not plotted those 
EMP stars known to be in binaries in Figure 1, but it is 
possible that many of these stars have as yet undetected 
companions. A uniform sample coverage of the CNO ele- 
ments would constrain this scenario better. 

In order to assess the implications of the nucleosynthetic 
data in EMP stars for reionization, we follow the approach 
in Madau & Shull (1996) and Venkatesan & Truran (2003). 
We calculate the ionizing efficiency, ??Lyc, which is essen- 
tially a conversion factor defined as the ratio of the energy 
produced in rest mass of metals, Mzc 2 , and the energy 
released in the H-ionizing continuum: 

VLyc = E(hv > 13.6eV)/(Af z c 2 ) (1) 

The first stars are expected to be metal-free in com- 
position, a factor that critically influences their structure 
and leads to the emission of significantly harder ionizing 
radiation relative to low-Z stars (Bromm et al. 2001b; 
Tumlinson et al. 2003; Schaerer 2002). Venkatesan & Tru- 
ran (2003) showed that r^yc is a strong function of stellar 
metallicity and increases by a factor of 10-20 from Z — Z & 
to Z = stars, owing to the boosted ionizing flux and 
substantially reduced metal yield from massive stars with 
decreasing stellar metallicity. The latter trend is a result 
of the more compact structure of metal-poor stars, which 
directly leads to greater element fallback onto the remnant 
after the SN explosion. 

The stellar IMF is assumed to have the form, <fi(M) = 
cf)oM~ a , which is normalized over the mass range as, 
J dMM(j)(M) = 1. We consider two primordial IMFs, 
one that has 100-1000 M Q stars with a flat IMF slope 
(hereafter very massive stars; VMSs), as some theoretical 
studies indicate (Omukai & Palla 2003; Bromm et al. 2002; 
Nakamura & Umemura 2002), and a present-day IMF with 
1-100 Mq stars where the IMF slope has a Salpeter value, 
a = 2.35. We use the ionizing spectra from Bromm et al. 
(2001b) and Schaerer (2002) for VMSs, and the results 
from Tumlinson et al. (2003) for the 1-100 M IMF. We 
take the metal yields from Woosley & Weaver (1995) for 
10-40 M Q stars, which end their lives as Type II SNe, and 
from Heger & Woosley (2002) for stars with masses 140- 
260 Mq which explode entirely as pair-instability SNe. All 
other massive stars are assumed to collapse into black holes 
with no contribution to the nucleosynthetic output. Vary- 
ing the upper and lower mass limits in the VMS IMF will 
not significantly impact our conclusions, owing to the nar- 
row mass range of element synthesis from VMS SNe, and 
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to the near-identical ionizing photon production per stellar 
baryon in all stars with masses > 300 Mq (Bromm et al. 
2001b). 

We do not include cases where the conversion efficiency 
is derived including the metal yield from intermediate- 
mass stars, as i]^ yc remains approximately the same for 
oxygen and exactly the same for Si. For carbon, rfi, yc de- 
creases by a factor of a few with intermediate-mass stars 
owing to their large yield in C. For the purposes of this 
paper, we wish to track the growth of ionizing radiation 
which is dominated by the massive stars in the IMF. We 
do not include the yield from Type la SNe in our calcula- 
tions. Recent studies indicate that such SNe may not occur 
in metal-free stellar populations (Kobayashi et al. 1998), 
although this is a sensitive function of the SN mechanism 
itself and of the unknown binarity in Z = stars. 

As shown in Venkatesan & Truran (2003) and Miralda- 
Escude & Rees (1997), calculations of riLyc can be used to 
compute the number of ionizing photons per baryon gen- 
erated in association with a metallicity that is observed 
in a given system. These papers focussed on the IGM, 
and pointed out that the universe first generates about 10 
ionizing photons per IGM baryon for the observed IGM 
metallicity of ~ 10~ 25 Z Q . This estimate is however sub- 
ject to the uncertainties of star formation efficiency and 
the escape fraction of ionizing radiation from early galax- 
ies. Here, we examine the ionizing efficiency in the en- 
vironment of EMPs, and calculate N 7 /Nb, the number 
of ionizing photos created per baryon in stars, which is 
not related to the above highly uncertain astrophysical 
parameters. Following Venkatesan & Truran (2003), for 
an element i, N 7ti /N b = i% yC}i x (1 GeV / {E Lyc )); 
we assume the IMF-averaged energy of a Ly-continuum 
photon, (E-Ly C } — 27 eV (30 eV) for metal-free stars in a 
present-day (VMS) IMF. 

3. RESULTS AND DISCUSSION 

For our calculations here, we take the values of 77Lyc 
for carbon and oxygen computed in Venkatesan & Truran 
(2003), and compute them independently for silicon. We 
find that for the 1-100 Mq Salpeter-slope IMF, rj^ yc has 
values of 0.48, 0.2 and 1.34 for C, O and Si respectively; 
for a VMS IMF, the corresponding values are 0.098, 0.01, 
and 0.057. Using these numbers, we show the results of 
our calculations in Figure 2, where the number of ioniz- 
ing photons created per baryon in EMP stars, Nj/N^, is 
shown as a function of [Fe/H] for the elements C, O and 
Si in three separate panels. In each case, the ionizing pho- 
ton contribution from a present-day and a VMS IMF are 
plotted. 

A number of trends seen in Figure 2 are of relevance to 
the first-stars field. First, VMS make a low but persistent 
contribution to the overall ionizing photon budget in asso- 
ciation with the metals locked in EMP stars, typically an 
order of magnitude less than the 1-100 Mq IMF. This is 
because VMSs are in general less efficient at generating to- 
tal ionizing radiation per unit metal yield relative to other 
metal-free stellar populations, owing to the large metal 
production from VMSs. Second, the slow buildup of met- 
als and the associated ionizing radiation is apparent, for 
each element and for each IMF. This trend appears most 
obvious for oxygen, although there is a significant scatter 



in the correlation between iV 7 /iVb and [Fe/H] at any given 
[Fe/H]. As we stated earlier, the notable exceptions are 
the two stars at the lowest values of [Fe/H]. 

Third, we show in each panel the approximate criterion 
to generate an optical depth in the cosmic microwave back- 
ground (CMB) of about 0.1, consistent with data from the 
Wilkinson Microwave Anisotropy Probe (WMAP; Spergel 
et al. 2006), of about N 7 /N h ~ 34,000 (Tumlinson 2006; 
Tumlinson et al. 2004). This number is derived from inter- 
polating between cases studied by Tumlinson et al. (2004) 
where the lifetime-integrated ionizing photon production 
from various stellar populations and IMFs were calculated 
and used as inputs in detailed cosmological reionization 
models. Interestingly, the EMP data indicates that this 
line is crossed at least once and perhaps twice between -4 < 
[Fe/H] < -3 by a present-day IMF. This is best seen in the 
panel corresponding to Si, where a clear dip in Nj/N^ oc- 
curs between -3.6 < [Fe/H] < -3. The VMS IMF appears 
to never meet this criterion over the metallicity ranges 
considered here, consistent with the role of VMSs in IGM 
reionization predicted by Venkatesan & Truran (2003). 
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Fig. 2. — The minimum number of ionizing photons per baryon 
in stars that is created in association with specific metals observed 
in MP halo stars. Top, middle and bottom panels correspond to the 
elements C, O, and Si. The thick dashed black line in each panel cor- 
responds to the value required to roughly match the WMAP data 
for an optical depth of 10% in the CMB. All cases are for metal-free 
stars; each panel displays values for the 1-100 Mq (diamond points) 
and VMS (triangle points) IMFs considered in this work. See text 
for discussion. 

Conducting a metal census of the universe from EMP 
stars and from the IGM is important in both environ- 
ments, despite their differing systematics. The IGM con- 
tains nearly all of the baryons at high redshifts but the 
former case permits a detailed star-by-star analysis of the 
cosmic enrichment history. By working with nucleosyn- 
thctic data in EMP stars rather than the IGM, we have not 
assigned a global fixed metallicity with the appropriately 
scaled individual element values according to their solar ra- 
tios. Rather, we have derived parallel constraints for reion- 
ization from a number of elements independently, which is 
a considerably stronger approach. In addition, observa- 
tions of EMP stars relative to IGM absorbers provide sig- 
nificantly larger element coverage in the well-constrained 
environment of stellar atmospheres, a situation that will 
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greatly improve in the near future. 

Although we naively apply the same factor for ionization 
efficiency for [Fe/H] > -3, r/Lyc is likely to be closer to the 
values associated with metal-poor stars, typically a factor 
of a few to 10 lower, owing to their increased metal yield. 
Depending on the individual metal abundance, this could 
lead to a less steep rise or flattening in the behavior of 
N 7 /Nb at [Fe/H] > -3. This prediction is somewhat more 
speculative, as the integrated contribution from many sub- 
sequent stellar generations of varying metallicities will be 
difficult to distinguish in this environment. 

The notable exceptions to most of the trends noted 
above are the two most Fe-poor stars at [Fe/H] ~ -5.5. At 
such low metallicities, one would expect that these stars 
reflect enrichment from a single SN (Tumlinson 2006). 
These stars are, however, highly over-enriched in several 
elements including C, O and Si, to the degree that they 
are actually not metal-poor, but rather at metallicities of 
nearly a tenth solar. How these stars came to acquire 
such skewed element abundances relative to solar ratios 
remains puzzling at present. Some possibilities include 
pollution from a binary companion, from a subluminous 
hypernova with fine-tuned ratios of mixing and fallback 
(Iwamoto et al. 2005), or the effects of stellar rotation and 
ensuing mass loss which could dramatically boost the pro- 
duction of CNO elements (Chiappini et al. 2006). The 
ratio of these elements to iron or relative to each other 
is sometimes used as a chronometer in other astrophysi- 
cal environments. For the most Fe-poor stars, however, it 
is clear that ratios such as [Fe/O] or [C/O] are not eas- 
ily related to time. These should normally increase with 
age of a stellar population as the stars of progressively 
lower mass evolve off the main sequence, but as we saw 
in Figure 1 for EMP stars, [C/O] achieves its highest val- 
ues at the lowest [Fe/H]. Perhaps these stars were created 
in atypical environments as the IGM approached reioniza- 
tion, as seen from the generation of ionizing radiation in 
Figure 2 - the most Fe-poor stars appear to have formed 
at epochs when, for a metal-free present-day IMF, Nj/Nt, 
nearly equals or exceeds the criterion from WMAP, at 
least for C. One interesting explanation of this, in addition 
to those above, could be that the most Fe-poor halo stars 
are also true second-generation stars: ones that formed 
through the dust grain cooling pathway (Z tI ~ 1O~ 5 Z ) 
rather than through metal line cooling using C, O or Si 
(Ztr ~ 10 Zq). The unusual element trends in these 
stars could then be explained through the element seg- 
regation predicted in a scenario involving dust transport 
in early galaxies. Here, conditions in primordial SN rem- 
nants result in the radiative transport of dust in the hot 
radiation field from metal- free stars, leading to selective 
element enhancement from compounds such as graphites 
and silicates in the cooling gas in SN shells (Venkatesan 
et al. 2006), which are likely the sites of EMP star forma- 
tion. 

What do the results in this paper imply for cosmology? 
Those EMP stars which have observational age estimates 



are about 12-13 Gyr old (Sneden et al. 2003; Beers & 
Christlieb 2005). In the context of our results, it appears 
that most of the EMP stars formed close to the end of 
or soon after cosmological reionization. That reionization 
occurs in the first Gyr of the universe (by z ~ 6) is no sur- 
prise, as this is already clear from Gunn-Peterson studies 
of the high-z IGM (Becker et al. 2001) and from current 
CMB data. A more interesting result of translating the 
EMP star abundances in C, O and Si into the accompany- 
ing production of ionizing photons per baryon in stars is 
the clear buildup of ionizing radiation until reionization is 
accomplished at [Fe/H] - -3, [Si/H] ~ -3.5 and [O/H] ~ 
-2.5 ([C/H] ranges from -4 to -2 at these Fe- metallicities). 
At this point, the initial stellar population (presumably 
metal-free) must turn off, as seen from the altered ele- 
ment abundance patterns at higher [Fe/H], leading to a 
new pattern of nucleosynthesis from different SNe. Cayrel 
et al. (2004) argue for the presence of a primordial SN 
population up to [Fc/H] or [Mg/H] values of about -3, 
with a combination of primordial and nonprimordial SNe 
at higher metallicities. Such a truncation in the mode 
of star formation is most likely attributed to the critical 
metallicity being reached in the gas in which EMP stars 
formed, signalling an end to the first-stars epoch. Given 
certain conditions of density and temperature, Z tI can be 
higher in iron than in carbon, oxygen, or silicon. This 
is not consistent, however, with the overall trend of en- 
hanced C, O and Si at low [Fe/H]. A naive projection of 
the above metallicities in Fe, Si, O and C at which reion- 
ization "occurs" in Figure 2 implies that EMP stars should 
have formed in gas of number densities > 10 3 cm -3 (San- 
toro & Shull 2006). 

A more puzzling coincidence is why the end of the first- 
stars phase should occur at about the epochs correspond- 
ing to cosmological reionization. This points towards the 
possibility of strong feedback between intrahalo metal rein- 
corporation timescales and halo-IGM transport of ionizing 
radiation. We return to this issue in a subsequent work, 
where we will quantitatively relate the reionization epoch 
to the transition metallicity in early galaxies. Clearly, 
EMP stars provide an important window to study the 
transition from first- to second-generation star formation 
and the physical conditions at which this occurs in pri- 
mordial galaxies. The novel result derived here is that the 
abundances in such stars show a strong relation to cosmo- 
logical reionization. Targetted observations are required 
with more uniform element coverage of EMP stars, partic- 
ularly the C-rich stars, in the critical range of -6 < [Fe/H] 
< -4. This range currently has very little data, and will 
be useful in closing the loop in our understanding of these 
important cosmic events at high redshifts. 
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